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' I i: A 1+3 dimensional solution of relativistic hydrodynamics is analyzed in this 

^1 paper. Momentum distribution and other observables are calculated from the 

solution and compared to hadronic measurements from the Relativistic Heavy 
Ion Collider (RHIC). The solution is compatible with the data, but only the 
freeze-out point of the evolution is determined. Many equation of states and 

I i initial states (initial temperatures) are valid with the same freeze-out distri- 

bution, thus the same hadronic observables. The observable that would dis- 
tinguish between these initial temperatures is momentum distribution of pho- 
^ tons, as photons arc created throughout the evolution of the fireball created 

in RHIC collisions. The PHENIX experiment at RHIC measures such data 

00 via low invariant mass e+e^ pairs. Average temperature from this data is 

T = 221 ± 23(stat)±18(sys) MeV, while a model calculation with initial tem- 
perature Tj„jt = 370 MeV agree with the data. 

— «^ Keywords: heavy ion collisions, hydrodynamics, equation of state, temperature, 

hadron spectra, photon spectra 

^ — I 

^ 1. Perfect fluid hydrodynamics 

In the last several years it has been revealed that the strongly interacting 
^ Quark Gluon Plasma produced^ in the collisions of the Relativistic Heavy 

Ion Collider (RHIC) is a nearly perfect fluid,^ i.e. it can be described with 
perfect fluid hydrodynamics. 

Perfect fluid hydrodynamics is based on local conservation of entropy 
or number density (n), energy-momentum density (T''^). The fluid is per- 
fect if the energy-momentum tensor is diagonal in the local rest frame, i.e. 
viscosity and heat conduction are negligible. This can be assured if T'"' is 
chosen as T^'^ = {e+p)u^u" —pg'^", where u** is the flow field in the fluid, e 
is energy density, p is pressure and g^'' is the metric tensor, diag(l,-l,-l,-l). 
The conservation equations are closed by the equation of state, which gives 
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the relationship between e and p. Typically e = np is chosen, where the 
proportionality "constant" k may depend on temperature T, which in turn 
is connected to the density n and pressure p via p = nT. 

The exact, analytic result for hydrodynamic solutions is, that the 
hadronic obscrvablcs do not depend on the initial state or the dynami- 
cal equations separately, just through the final state. "^'^ Thus if we fix the 
final state from the data, the equation of state can be anj^hing that is 
compatible with the particular solution. This is the framework of several 
hydro solutions as detailed in the next paragraph. 

Many solve the above equations numerically, but there are only a few 
exact solutions. Historically the first is the implicit 1+1 dimensional ac- 
celerating solution of Landau and Khalatnikov.^^^ Another renowned 1-1-1 
dimensional solution of relativistic hydrodynamics was found by Hwa and 
Bjorken:^"^*^ it is simple, explicit and exact, but accelerationless. 

Important arc solutions^^'^^ which are explicit and describe a relativistic 
acceleration, i.e. combine the properties of the Landau-Khalatnikow and 
the Hwa-Bjorken solutions. With these one can have an advanced estimate 
on the energy density, but investigation of transverse dynamics is not 
possible by these solutions. 

The only exact 1-1-3 dimensional relativistic solution, from which ob- 
servables like momentum distribution, correlation function and elliptic flow 
were calculated'* is the one in ref.^"* Observables from this solution were 
computed and compared to data in ref.^ 



2. The analyzed solution 

The analyzed solution^* describes an ellipsoidally symmetric expansion. 
The ellipsoids are given by constant values of the scale variable s: 

x{tY Y{ty z{ty' ^ ' 

here X{t), Y{t), and Z(t) are time dependent scale parameters (axes of the 
s = 1 ellipsoid), only depending on the time t. Spatial coordinates are r^, 
ry, and r^. The velocity-field is described by a Hubble- type expansion: 

( . X(t) Y{t) Z(t) \ 

where x means the four-vector {t,rx,ry,rz), and X(t) = dX(t)/dt, similarly 
for Y and Z. The X{t) = Xq, Y{t) = Yq, Z{t) = Zq (i.e. all are constant) 
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criteria must be fulfilled, ic. the solution is accelerationless. This is one of 
the drawbacks of this solution. 

The temperature T{x) and number density n{x) are: 



where t is the proper time, s is the above scaling variable, i>(s) is an arbi- 
trary function, while no = n|s=o,T=To, = T\s=o,t=to and po = p\s=o,t=to 
with Po = noTo (hence p does not depend on the spatial coordinates only 
r). Furthermore, To is the time of the freeze-out, thus Tq is the central 
freeze-out temperature. The parameter k is arbitrary, i.e. any value of k 
yields a solution. The function h'(s) is chosen as: 



where b is then the temperature gradient. If the fireball is the hottest in 
the center, then 6 < 0. An example time evolution of the temperature 
distribution is shown in fig. 1 

3. Hadronic observables 

The picture widely used in hydro models is that the pre freeze-out (FO) 
medium is described by hydrodynamics, and the post FO medium is that 
of observed hadrons. In our framework we assume that the freeze-out can 
happen at any proper time, e.g in case of a self-quenching effect or if the 
phase space evolution is that of a collisionless gas. See details in ref.'' The 
hadronic observables can be extracted from the solution via the phase-space 
distribution at the FO. This will correspond to the hadronic final state or 
source distribution S{x,p). 

We do not need to fix a special equation of state, because the same final 
state can be achieved with difi'erent equations of state or initial conditions."^ 
Thus the hadronic observables do not restrict the value of k. 

The first calculated observable is invariant transverse momentum dis- 
tribution Ni{pf) of a particle with mass to:"* 




(3) 




(5) 



K5) = e-''^/2, 



(6) 
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Fig. 1. Temperature distribution in the transverse plane (x-y) is shown for various times 
with an example parameter set. The fireball is the hottest in the center and it cools down 
as time elapses. 



with the following auxiliary quantities: 



iV=AAno('^V''. (7) 



nit 



> (9) 



Furthermore, ,Ty, are the effective temperatures, i.e. inverse logarith- 
mic slopes of the distribution: 

T^=T.+^;:f^y (10) 

o(To - mt) 

T.=To+-;^?^, (11) 

6(To - mt) 

T.=To + ^*^, (12) 
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where Xq, Yq and Zq are the (constant) expansion rates of the fireball, Tq 
its central temperature at FO and b the temperature gradient. 

We also calculate the elliptic flow, which describes the azimuthal asym- 
metry of the momentum distribution:^ 

where Iq, and Ii are the modified Bcsscl functions while 

»=#f^-^)- M 



4mt \Ty 

See details of the calculation in ref.** The formula for V2 gives back pre- 
viously found formulas of non-relativistic solutions^^ and relativistic solu- 
tions. ^^'^'^ Also the formula for Ni{pt) is similar to results of the previously 
mentioned papers. 

Third observable we calculate is the two-particle Bose-Einstein (HBT) 
correlation radii of identical bosons: 

2 _ TpT^jT^ - To) 

' M^T, ' ^^^^ 

2 _ ToT§{Ty - To) 
- ' (Id) 



MtT, 




ToT§{Ty- 


To) 


MtTy 






To) 



where Mj is the transverse mass belonging to the average momentum K = 
0.5(pi +P2) of the pair, which is (at mid-rapidity) Mt = 0.5 (to(,i -|- mt^2 )■ 
The rnt^i, and mt,2 quantities are the transverse masses, the T^, Ty, and 
Tz are the effective temperatures belonging to the average momentum (i.e. 
here = T^ImJ- The calculations are detailed in ref.^ To compare the 
HBT radii with the data the Bertsch-Pratt^^ frame is to be used. It has 
three axes: the out is the direction of the average transverse momentum of 
the pair, the long direction is equal to the direction z, and the side direction 
is orthogonal to both of them. The result for i?out, -Rside and i?iong is: 

^2 _|_ j^2 

-^out = -^side = (18) 

Aong = Rl- (19) 

Clearly in this solution the out and side radii are equal. This can be at- 
tributed to the instantaneous freeze-out; a non-zero freeze-out duration 
would make R^^^ bigger by a term of At^p^/E"^. Supported by the data. 
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we use the At = approximation in our solution, which corresponds to 
instantaneous freeze-out. 

An important consequence of the above results is that neither spec- 
tra nor elliptic flow nor correlation radii depend on the EoS itself, only 
through the final state parameters. If we determine for example To, the 
freeze-out central (at the center means here Tx — Vy — rz = 0) tem- 
perature, K or the initial temperature Tjnitiai still cannot be calculated. 
We only know that they are connected through the consistency condition 
^initial = To(To/Tinitiai)^^'', scc cq. (4), i.e. they can be co- varied (softer 
EoS requires smaller initial temperature for a given freeze-out proper-time). 
Thus K or Tinitiai has to be determined from another measurement, e.g. the 
spectrum of thermal photons. 

4. Comparing the hadronic observables to RHIC data 

The above results were compared in ref.^ to PHENIX data of 200 GeV 
Au-|-Au collisions. Above formulas were fitted to describe spectra and HBT 
positive pion data^°'^^ (0-30% centrality) and elliptic flow data^^ for tt^, 
K"^, p and p particles (0-92% centrality). 

The fit results are shown in fig. 2, see details of the fit in ref.^ Important 
are the following: central freeze-out temperature Tq is around 200 MeV for 
both datasets (with an error of 7 MeV), and the fireball is colder away 
from the center. The expansion eccentricity is positive, it tells us that the 
expansion is faster in-plane. Because of the Hubble-flow this means that the 
source is in-plane elongated, similarly to the result of ref.^^ The freeze-out 
happens at a proper-time of tq = 7.7 ± 0.8 fm/c. 

Our fit parameters describe the fireball at the freeze-out. However, the 
solution is time-dependent, most importantly the temperature depends on 
time as described by eq. (4) . We plotted the time-dependence of the central 
temperature in fig. 3 for several values of k, i.e. several EoS'. From this, 
assuming for example an average k of 10, see ref.,^^ one can also calculate 
the initial central temperature of the fireball based on eq. (4): 



This yields 370 MeV at tinitiai=l fm/c (note that t = r at the center). 

One still would like to determine the value for k and the initial tem- 
perature from experimental observables. The key are thermal photon in- 
variant momentum distributions, because photons are not suppressed by 
the medium, they are produced according to the local temperature at all 




(20) 
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times. If one could measure the yield of thermal (direct) photons, one could 
compare it to time-integrated hydrodynamic results and determine EoS and 
initial temperature. In the next section thus experimental results on direct 
photons are reviewed. 




0.18 




200 400 600 800 1000 1200 1400 1600 1800 



p, [MeV/c] 

Fig. 2. Fits to 0-30% centrality PHENIX Au+Au spectra^" (top left) HBT radii^^ (top 
right) and 0-92% centrality PHENIX Au-(-Au elliptic flow^^ (bottom). See details of the 
fit in ref.* 



5. Direct photons at the PHENIX experiment at RHIC 

Thermal photons from the partonic phase are predicted to be the domi- 
nant source of direct photons for 1 < < 3 GeV/c in Au+Au collisions at 
the Relativistic Heavy Ion Collider (RHIC).^^ The measurement of direct 
photons in this kinematic domain is however very difficult due to the back- 
ground from hadronic decay photons. At PHENIX, an alternative approach 
is used.^^ The idea is that any source of high energy photons can also emit 
virtual photons which then convert to e+e~ pairs. 
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Fig. 3. Time dependence of the central temperature of the fireball, from eq. (4) is shown 
for different k values. In reality k may change with time, we show here the curves only 
for fixed k values. Assuming an average of k = 10^* one gets an initial temperature of 
370 MeV at tinitial=l fm/c, in agreement with PHENIX measurements.^^ 

The relation between photon production and the associated e+e" pair 
production is then^^^^^ 



Here a is the fine structure constant, me and m are the masses ol the 
electron and the e'^e'^ pair respectively, and 5 is a process dependent factor. 

Invariant mass distribution of e^e~ pairs was measured thus in 
PHENIX25.27 foj. < MeV/c^ and for 1 < < 5 GeV/c in Au+Au 
and p -\- p collisions at y/s^N = 200 GeV. Top panel of fig. 4 shows the 
measured^^ mass spectra of e+e^ pairs in p + p and Au+Au collisions for 
different ranges of pt, comparing them to expected yields from dielcctron 
decays of hadrons, calculated using a Monte Carlo hadron decay generator 
based on meson production as measured by PHENIX. The Au+Au data 
show a relatively high excess above the hadronic background, which indi- 
cates an the production of virtual photons in Au+Au collisions. PHENIX 
assumed that the excess is entirely due to internal conversion of direct pho- 
tons and deduce the real direct photon yield from the e^e~ pair yield using 



In order to calculate direct photon yields the /(m) = (1 — r)fc{m) + 
ffdiri'iTi) fit function was used to the mass distribution in each pt bin sep- 
arately, where fc{m) is the mass distribution from the hadronic decays 




(21) 



eq. (21). 



28 
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Fig. 4. The e+e~ pair invariant mass distributions in p + p and minimum bias Au+Au 
collisions from ref.^^ is shown in the top panel. Mass distribution of electron pairs mea- 
sured in Au+Au minimum bias events'^^ is shown for 1.0 < pt < 1.5 GeV/c in the 
bottom panel. The fit is explained in the text and originally in ref.'^^ 

(shown in the top panel of fig. 4 with sohd Unes), and fdiri'm) is the ex- 
pected shape of the direct photon internal conversion, and r is the fit pa- 
rameter. Such a fit is shown in the bottom panel of fig. 4 for one pt bin. The 
direct photon yield can be calculated then using the direct photon fraction 
r as dN'^'''"'\pt) = r X diV'"'^^"'^'™(pt) where diV'"'=i"'^'™(pf ) is the yield of 
all photons (direct and decay inclusive). The inclusive photon yield can be 

calculated as dArin<=l"«ivc ^ ^data ^ (^^cocktalY^cocktall)^ ^jjgj.g ^data 
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^cocktail g^j.g measured and cocktail e+e" pair yields and c^^cocktaii 
the yield of photons from the cocktail. 




Fig. 5. Invariant cross section (p + p) and invariant yield (Au+Au) of direct photons 
is shown, taken from ref.^^ The red points are from the analysis presented in ref.^^ and 
blue points are from.^®'^" Dashed black curves are modified power-law fits to the p + p 
data,^^ while the solid black curves are exponential plus modified power- law fits.25 The 
red dashed curve is the direct photon spectrum in central Au-|-Au collisions in ref.'^® 

The obtained direct photon spectra^^ are shown in fig. 5 and compared 
to similar data of refs.^^''^° The direct photon yields are consistent with 
a NLO pQCD calculation in p + p. The shape of the direct photon spec- 
tra above a binary collisions scaled p + p spectrum is exponential in pt , 
with an inverse slope T = 221 ± 23(stat)±18(sys) MeV in central Au-|-Au 
collisions. This temperature can be regarded as a time-average of the 
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temperature of the fireball, represents thus an experimental lower limit to 
the initial temperature. The shape of the thermal photon spectrum from 
a model calculation with initial temperature Tinu = 370 MeV agrees with 
the data.^^ 

6. Summary 

Exact parametric solutions of perfect hydrodynamics were long searched for 
in order to describe the matter produced in heavy ion collisions at RHIC. 
We extracted hadronic observables from the relativistic, 1+3 dimensional, 
ellipsoidally symmetric, exact solution of ref."'^^ We calculated momentum 
distribution, elliptic flow and Bose-Einstein correlation radii from the so- 
lution. We compared the results to 200 GeV Au+Au PHENIX data.^^-^^ 
The solution is compatible with the data. If using an experimentally de- 
termined average EoS of k « 10, our results yield approximately 370 
MeV at rinitiai=l fni/c, in agreement with recent PHENIX photon mea- 
surements.^^ The time-average of the temperature in these collisions is 
r = 221 ± 23(stat)±18(sys) MeV,^^ which is a lower limit of the initial 
temperature. From detailed comparisons with hydro dynamical models of 
direct photon emission, PHENIX concluded that Tinu = 300 MeV is the 
lowest possible initial temperature at 1 fm/c, that is consistent with such 
an average slope parameter. In the future we will compare hadron and 
photon observables simultaneously to the analyzed analytic hydrodynamic 
solution. 
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